BA-TH/03-463 
CERN-TH/2003-128 
hep-th/0306113 



Perturbations in a non-singular bouncing Universe 



o 

(N 
C 

Hi 

(N 



CO 



M. Gasperini 1 , M. Giovannini 2 and G. Veneziano 2 



1 Dipartimento di Fisica, Universita di Bah, 
Via G. Amendola 173, 70126 Bah, Italy 
and 

Istituto Nazionale di Fisica Nucleare, Sezione di Bah 
Via G. Amendola 173, 70126 Bah, Italy 

CO ' 

o 

' 2 Theoretical Physics Division, CERN, 

CH-1211 Geneva 23, Switzerland 

> 

X: 

^ , Abstract 

We complement the low-energy gravi-dilaton effective action of string theory with 
a non-local, general-covariant dilaton potential, and obtain homogeneous solutions de- 
scribing a non-singular (bouncing-curvature) cosmology. We then compute, both an- 
alytically and numerically, the spectrum of amplified scalar and tensor perturbations, 
and draw some general lessons on how to extract observable consequences from pre-big 
bang and ekpyrotic scenarios. 



In recent years string/M-theory has inspired new cosmological scenarios in which a long 
period of accelerated (growing-curvature) evolution, emerging from an almost trivial initial 
geometry, turns into a standard (decreasing-curvature) FRW-type cosmology, after going 
smoothly through a big bang- like event. There are by now several variations on this pre-big 
bang theme. Besides the original pre-big bang (PBB) scenario [1, 2, 3], based on the duality 
symmetries of string cosmology, new models incorporating brane and M-theory ideas have 
been proposed under the generic name of ekpyrotic (EKP) scenarios [4]. While different 
proposals differ in the way the scale factor behaves during the growing-curvature phase, 
they all share the feature of describing a bounce in \H\ (the absolute value of the Hubble 
parameter) or, in more geometrical terms, in the space-time curvature. The latter always 
starts vanishingly small, grows to a maximum at the would-be big bang, and then decreases 
again in the FRW phase. A common theoretical challenge to all these models is that of 
being able to describe the transition between the two regimes. 

At a more phenomenological level, instead, the challenge is to compute, in a reliable 
way, the final spectrum of amplified quantum fluctuations to be compared with present 
data on CMB radiation and large-scale structure. In the PBB case it was admitted early 
on [5, 6] that adiabatic- curvature perturbations had too large a spectral index to be of 
any relevance at cosmologically interesting scales (while being possibly important for grav- 
itational waves searches [7, 8]). Isocurvature perturbations (related to the Kalb-Ramond 
two- form) can instead be produced with an interestingly flat spectrum [9], but have to be 
converted into adiabatic-curvature perturbations through the so-called curvaton mechanism 
[10] before they can provide a viable scenario for large scale-anisotropies [11]. Proponents 
of the ekpyrotic scenario, while agreeing with the PBB result of a steep spectrum of tensor 
perturbations, have also repeatedly claimed [4, 12] to obtain "naturally" an almost scale- 
invariant spectrum of adiabatic-curvature perturbations, very much as in ordinary models 
of slow-roll inflation. These claims have generated a lot of heated discussion (see for in- 
stance [13]), with many arguments given in favour or against the phenomenological viability 
of EKP scenarios in the absence of a curvaton's help. The reasons for the disagreement 
can be ultimately traced back to the fact that the curvature bounce (hereafter simply re- 
ferred to as the bounce) is put in by hand, rather than being derived from an underlying 
action. This leaves different authors to make different assumptions on how to smoothly 
connect perturbations across the bounce itself, and this often results in completely different 
physical predictions. 

Our aim is to present a class of models where a regular bounce is derived from the 
field equations of a general-covariant, albeit non-local, action, and where perturbations can 
be studied from beginning to end. We are thus able to test, at least within the model, 
the assumptions made by different groups about how different perturbations should or 
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should not behave across the bounce. In this letter we only consider a simple gravi-dilaton 
model and present the main results, working in the string frame and skipping their detailed 
derivation. In a longer, forthcoming paper [14] we shall give full details of the calculations in 
both the string and Einstein frames (with identical physical conclusions) , and generalize the 
class of models to include a fluid component. We first present the action, the background- 
field equations, and a class of smooth, bouncing solutions. We then discuss, successively, 
the amplification of tensor and scalar perturbations from an initial spectrum of vacuum 
quantum fluctuations, and compare analytical and numerical results. We finally summarize 
our results and draw some conclusions. 

We recall that the homogeneous cosmological solutions derived from the tree-level, low- 
energy string effective action exhibit in general a curvature singularity, disconnecting the 
pre-big bang branch from the post-big bang one [15]. Such a singularity is expected to be 
removed by higher-loops and higher-curvature a' corrections (see for instance [3]), but it is 
known [16] that non-singular solutions can be explicitly obtained already at low curvatures 
[2] in the presence of an appropriate non-local effective potential. An example of such a 
possibility can be illustrated by using a (d+ l)-dimensional, general covariant gravi-dilaton 
effective action, which, in the string frame, reads: 



l — f d d+1 x^\ e~* [R + (V^) 2 + V(jp)} (1) 



2X a s 

(metric conventions: H . . .). Here A s = a/o 7 = Mf 1 * s the string length scale, and the 

potential V(Tp(x)), a local function of Tp, is instead a non-local function (yet a scalar under 
general coordinate transformations) of the dilaton owing to the definitions: 

d d+i 



V = V(e-n, e-*W = / yj\^\ e-^yVWMy) H<p{x) - <p(y)) . (2) 



Note that plays the role of a "reduced" coupling constant from d + 1 to + 1 space-time 
dimensions. We may thus expect V{e^ ip ) to go like some inverse power of its argument 
as this becomes large (i.e. in perturbation theory). We shall discuss in [14] how such a 
non-local potential may be induced by loop corrections in (higher-dimensional) manifolds 
with compact spatial sections; here, we take this simply as a toy model that avoids the 
singularity, while staying all the time at low energy/curvature. For a background manifold 
isometric with respect to d spatial translations, it is known that action (1) leads to field 
equations that are covariant under scale-factor duality [1] and, in the presence of the two- 
form background B^, also under global 0(d,d) transformations [16]. 

Variation of the above action with respect to g^ v and ip, though somewhat unusual, is 
straightforward and leads to the following field equations (in units of 2Ag _1 = 1): 
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R + 2VV - (Vtf + V - f + e-^gp/. - e-'V/, = 0, 
where (with a prime denoting differentiation with respect to the argument) 

h = hJ dd+ly ^y)\ y\^ iy) ) s m*) - ^)) , 

/ 2 = j d d+l yj\^y)\ ^d^{y)d»p{y) 5' (<p(x) - <p(y)) , 

and we have also introduced the induced metric and Laplacian: 

dp<pd v <p ~ 2 



(4) 



l[iv — 9jJ,v 



vv = >v*v>. 



(5) 



(6) 



For a homogeneous, isotropic and spatially flat background, we can set <7oo — 1) 9ij — 
— a 2 (t)8ij, ip = <p(t), and obtain e - ^ = e - ^^, where we have absorbed into ip the dimen- 
sionless constant — ln(J d d y/\g), associated with the (finite) comoving spatial volume. In 
such a case, the time and space components of Eq.(3), and the dilaton equation (4), lead 
to a set of equations already studied in [16, 2, 5]: 



Tp - dH 2 - V = 0, 

• 2 9 dV 
2^-lp 2 -dH 2 + V- — = 0, 



H — HTp = 0, 



(7) 
(8) 



where the third (dilaton) equation follows from the first two, provided ip ^ 0. 

As noticed in [16], the above set of equations admits, quite generically, non-singular 
solutions. To recall this point and give an explicit analytical example we note that the two 
equations (7) can be reduced to quadratures, i.e. 



H = me"?, 



= / d A< 



d + X V(m\') 



-1/2 



(9) 



where m is an integration constant, and A = m~ 1 e~ lp . If the function appearing above 
between square brackets has a simple zero, we obtain a regular bouncing solution [16]. 
Furthermore, if VX 2 —> at large |A| (corresponding to a potential generated beyond two 
loops), the asymptotic solutions approach the usual vacuum solutions at large \t\. 

Consider, in particular, the class of potentials 

2-l/n 



V{mX) = A a- (mX) 



-In 



-d 



(10) 



parametrized by the dimensionless coefficients a and by the "loop-counting" parameter n. 
For a > and n > 0, Eq. (10) leads to the general exact solution 

l/2n 



H = me* 



m 



a 



1 + (ami) 



2n 



(11) 
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As \t\ — > oo, the Hubble parameter behaves like H\t\ ~ a( 1 ~ 2n )/ 2n , so that the minimal pre- 
big bang solutions dominated by the dilaton kinetic energy are only recovered for a^ 1 /™) = 
d (see for instance [3]). If this condition is not satisfied, the background may still be non- 
singular, but the dilaton potential cannot be neglected, even asymptotically. 

For the perturbation analysis of this paper it will be sufficient to use as a toy model the 
simple regular background associated to the four-loop potential ViTp) = —V^e^. In this 
case the general solution, in d = 3 spatial dimensions, can be written as: 



a{r) 



T + vV 2 + 1 



I/n/3 1 

Tp = --ln(l + T 2 ) + v? , r = t/t Q , (12) 



where (po is an integration constant and £q 1 = e ipo \/Vo. The solution is thus characterized 
by two relevant parameters, tg 1 and e^ , corresponding, up to some numerical factors, 
to the Hubble parameter and string coupling at the bounce, respectively (without loss of 
generality we have set a(0) = 1 with t = 0, the time at which the bounce occurs). We 
stress that the possibility of regular backgrounds is not limited to the class of potentials 
illustrated in Eq. (10), and that regular bouncing solutions can also be obtained by adding 
to the action (1) fluid matter sources, as already pointed out in [2] and illustrated in [14]. 

The evolution equation of tensor (transverse and traceless) metric perturbations can be 
obtained by perturbing to first order the (i, j) component of Eq. (3). For each polarization 
we obtain, in Fourier space, 

h k — Tphk + u) 2 hk = , u) = k/a. (13) 

We shall concentrate our attention on the case of regular solutions approaching asymptot- 
ically the well- known minimal gravi-dilaton model with negligible potential. The gravita- 
tional wave spectra arising in this case are characterized, in general, by a steep spectrum 
{tit > 1), up to logarithmic corrections. In fact the asymptotic solution of Eq. (13) for 
large wavelengths can be written as 

h k = A k + B k f ^rdt, (14) 

where t ex ~ uj~ 1 denotes the time at which the perturbation exits the horizon. Since in the 
case under consideration e v a^ 3 ~ for t — > — oo, we must expect a Inojt growth of h k , 
leading, ultimately, to In k corrections in the power spectrum. 

For an accurate camparison of analytical and numerical results, we shall restrict our 
attention to the specific regular bouncing solution (12). In this case, the evolution equation 
for the canonical normal mode \x k = ae~ v ^ 2 h k has asymptotic solutions with normalization 
to an initial vacuum fluctuation spectrum given in terms of Hankel functions of index zero. 
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The application of the standard matching procedure leads to a Bogoliubov coefficient (see 
below for its exact definition) given by 



\[3k\ 2 = e 1 + e 2 ln 2 (k 1 /k), (15) 

where ei and e 2 are numerical factors of order 1, and k\ is a characteristic momentum scale. 
The associated power spectrum is 

where the normalization H(Q)/Mp(Q) ~ (A s /io) e^^l 2 is controlled by the ratio between 
the curvature scale and the Planck mass at the bounce t = 0. 

For the numerical calculation, it is convenient to work with the rescaled variable h = 
e _lp / 2 a 3 / 2 ft, = s/afi. In order to compute the amplification, we impose for t — > — oo the 
quantum mechanical initial conditions 

Mfc = -£=e- ifc " h k = -Le^/^*, (17) 



where we recall that we are using units in which M 2 = 2. The action of the gravi-dilaton 
background will produce a mixing in the positive and negative frequency modes so that, for 
t — ► +oo, the solution can be parametrized as 

h = ^= (afce"*/"* + &e < / w *) . (18) 



2w 

In a second quantized approach |/3fe| 2 is nothing but the number of produced gravitons of 
given momentum k. The energy density of the produced gravitons will then be, up to 
numerical factors, / din A; k A \j3 k \ 2 . Thus, using the theoretical expectations given above, 
this corresponds to an energy spectrum going as A; 4 (up to logarithmic corrections). The 
numerical integration in terms shows that both h and h pass regularly through the bounce; 
this is illustrated in Fig. 1, where we report the tensor power spectrum, \Shk\, in the region 
around t = 0, clearly in qualitative agreement with Eq. (16). The absolute normalization of 
the tensor power spectrum has been fixed by imposing H(0)/Mp(0) = 10~ 2 . 

The mixing coefficients of Eq. (18) can be obtained in terms of the asymptotic (r — ► +oo) 
values of the real and imaginary parts of h, via the expressions: 



K| 2 + |/3 fc | 2 = (io\h k \ 2 + - 



hk — 



(19) 

a k \ 2 -\Pk\ 2 = i(h k k k -h k h k y (20) 
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Figure 1: Behaviour of the tensor power spectrum \Shk\, for different values of k = kto, in 
the region around the bounce. In the present and in the following figures the log denotes 
the logarithm in basis ten. 

We expect Eq. (20) to be identically equal to 1 thanks to the conservation of the Wronskian. 
On the other hand, Eq. (19) should approach a (^-dependent) constant only at late enough 
times. These two expectations are perfectly fulfilled, as illustrated in Fig. 2, which thus 
represents a highly non-trivial consistency check of our numerical procedure. 

Clearly, smaller k re-enter the horizon (and thus reach their asymptotic value) later. 
The asymptotic value gives, for each k, the sought after Bogoliubov coefficient. The results 
are plotted in Fig. 3 (crosses) and fitted to the theoretical expectation (15). The fit is very 
good and the parameters ei, 62, k\ are determined according to: 

|/3 fc | 2 = 0.46 ln 2 (fc/A;i) -4.84, k t = 6.68/V (21) 

We have not yet attempted to make a theoretical estimate of 62- 

Unlike tensor perturbations, scalar fluctuations are, in general, neither gauge- nor frame- 
independent. Gauge-invariant perturbations can be defined, as usual, in each frame. Among 
these, the spatial curvature perturbation on uniform dilaton hypersurfaces 1Z, and the 
Bardeen potentials \P and $ are particularly important for CMB phenomenology. Denot- 
ing by S(p = x the dilaton perturbation, and introducing scalar metric perturbations in 
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Figure 2: Asymptotic behaviours of the mixing coefficients given in Eqs. (19) and (20) for 
different values of k = kto. 

conformal time via the standard expressions 

5g 00 = 2a 2 cf>, bgij = 2a 2 (ipdij - didjE), Sg i0 = a 2 d { B , (22) 
we may express H, $ and <5 as: 

n = -i>-^- X , * = ip + H(E'-B), <& = </> — H{E' — B) - (E' — B)', (23) 

f 

where primes denote derivatives with respect to conformal time and Ti = a' /a. Introducing 
the Einstein frame, with scale factor cle = ae~' fi / 2 , we have 4>e = (f>s — \x-, = V's + \x 
(while B and E coincide in the two frames), and we can easily prove that 1Z and + are 
frame- independent, while (\& — $) s = (vf — &)e — \ + ip'E'. 

We shall work in the convenient uniform dilaton gauge x = 0, and also set B = 0, so 
that 1Z = —tp and ^e = tP + TCeE', where TLe = H — <p'/2. This gauge is particularly useful 
since all the perturbations variables in (22) coincide in the two frames. Another advantage 
is that, by perturbing the background equation (3), we find in this gauge 5^ = 0, 5e~ v = 0, 
5I\ = 0. The perturbations of the background equations thus simply lead to the following 
set of evolution equations: 

3(<p' - 2H)ip' + (6H<p' - if' 2 - 6H 2 )(j) + 2V 2 ip - (<p' - 2H)V 2 E' = 0, (24) 
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Figure 3: The numerically determined values of |/3fc| 2 are fitted to theoretical expectations 
for tensor (crosses) and scalar (stars) perturbations. 



(<ft - 2^)0 = 2V', 
E" + (2H - <p')E' + i) 



0, 



(25) 
(26) 



following, respectively, from the (00), (Oi) and (i ^ j) components of the perturbed Eq. (3). 
Note that the last equation, upon using He = H — <p'/2, is nothing but — &)e = 0, and 
that the additional equations obtained from the (i = j) component and from the dilaton 
equation (4) are redundant. 



We may combine the above system to get the interesting equation 
1/ = (v/ ~fV [y - 2H)E> - 2*] , 
which turns out to be exactly the standard evolution equation for 1Z: 



n' 



if 



(27) 



(28) 



Although this equation is formally the same as the standard equation of ordinary cosmo- 
logical perturbation theory, the relation between ip' 2 to (H E — H' E ) is different in our case. 
Indeed, it can be easily checked, by transforming the background equations (7), (8) to the 
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Einstein frame, that: 

4(H' E -H 2 E ) + <p' 2 + ^a 2 e e* = 0. (29) 

This shows that the usual identity ip' 2 oc (Ti' — H 2 ) is not satisfied in our case. And indeed, 
in our class of bouncing solutions <p' , unlike 7i' E — Ti. E , never vanishes. As a consequence, 
the pre-factor on the right-hand side of Eq. (28) does not diverge during the bounce; we can 
thus infer, already at this stage, that 7Z is likely to be conserved on super-horizon scales. 

Let us now turn to the specific study of the background model already analyzed in the 
case of the tensor perturbations. On the basis of simple estimates based on the asymp- 
totic behaviour of the solutions prior to the bounce and after, we expect the amplification 
coefficients for the scalar and the tensor modes of the geometry to be very similar, both qual- 
itatively and quantitatively. Indeed, after some lengthy but trivial algebra, the following 
decoupled evolution equation for 1Z can be obtained from Eqs.(24)-(26): 

1Z" k + 2^TZ' k + k 2 C 2 lZ k = 0, (30) 

where 

On the other hand, in order to enforce the correct quantum normalization, we note that 7Z 
is related to the canonical normal mode of scalar perturbations, v k = z7Z k . For large \t\, 
V — > 0, z oc ae^^l 2 and we are led exactly to the same equations as for tensor perturbations, 
the same asymptotic solutions, and thus the same spectrum for SR k . 

These expectations will now be checked numerically. We first notice that it is impossible 
to follow the evolution of perturbations directly through the canonical variable v (or the 
other often used variable u), since these become singular at the (Einstein- frame) bounce. 
This is not a problem, however, since v is only needed at very early or very late times, for 
normalization purposes. We can instead integrate directly a first-order system of differential 
equations, involving ip = — TZ and E' , which is completely regular throughout. Going over 
to (string- frame) cosmic time, and using the constraint (25) into Eqs. (24) and (26), we 
can eliminate <fi. We then obtain: 

ip k = A k {t)i> k + B k (t)£ k , 4 = C k (t)ip k + D k {t)S k , (32) 

where £ k = a 2 E k , and 

2(0 -2H) 2 £> f_i_\ {<P-2H\ 2 2 



A k {t) = ^— B k (t) = - 



Auj 2 

C *( t ) = ^r- 1 > D k (t) = <p-H -A k (t). (33) 
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Figure 4: Behaviour of the power spectrum of 1Z, for different values of k = kto. 

This system can be solved by imposing quantum mechanical initial conditions for the fluc- 
tuations of the curvature. Using v k = zlZ k we expect, asymptotically 

U k ~ ±-±={a k e- l f" dt + f3 k e l f" dt ), (34) 

where w = y/az. 

Following the same lines the discussion as were developed in the case of the tensor modes 
of the geometry, we can obtain, in the asymptotic region after the bounce, the appropriate 
combinations that determine the mixing coefficients: 

\a k \ 2 -\Pk\ 2 = iv?(n%n k -n k n%), 

\a k \ 2 + \(3 k \ 2 = w 2 L\K k \ 2 + ^\YTZ k + 7t fc | 2 ) , (35) 

where Y = z/z. The numerically computed power spectrum of TZ is reported in Fig. 4. 
In full analogy with the procedure discussed before, we have also plotted the quantities 
appearing in Eqs. (35), checking that the first one gives identically 1, and that the second 
approaches a limiting value at late times. We can thus extract the corresponding value of 
|/3fc| 2 and fit it to the theoretical expectation. The results of this analysis are reported in 
Fig. 3, together with the already discussed results for the tensors. The fit for the scalar 
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Figure 5: Time evolution of the power spectrum of the (Einstein-frame) Bardeen potential 
(thin and dashed curves). The bold curve represents the behaviour of an (appropriately 
rescaled) 1Z mode. 

case (starred points) gives: 

\(3 k \ 2 = 0.46 \n 2 {k/k 2 ) - 2.22, k 2 = 2.2/t . (36) 

This formula is very similar to the one obtained in the case of tensors perturbations. 

We now turn our attention to the behaviour of the (Einstein-frame) Bardeen potential. 
Its time evolution is shown in Fig. 5, where it is clearly visible that becomes much larger 
than one for small kt], signalling a breakdown of perturbation theory in the longitudinal 
gauge (a point already emphasized by the first paper quoted in Ref. [13]). The two Bardeen 
potentials in the string frame, though unequal, exhibit similar pathological behaviours. We 
also see that, although all variables are regular across the bounce, the values of \I/ and 
unlike those of 1Z, change drastically across the bounce. Thus, assuming continuity of ^> 
and of its derivative in simplistic models for the bounce is very dangerous. 

Let us finally turn to another popular variable, the spatial-curvature perturbation on 
constant-density hyper surf aces, £. It can be shown from its definition that it is not frame- 
independent, i.e. that 

Cs = -i> B - ^jSp s + (e = ~1pE - ^T^PE- (37) 
Ps Pe 
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In any case, in our gauge, Qe and Q s differ by 1Z by terms proportional to (p. The latter, 
however, is very suppressed at large scales (cf. Eqs. (25) and (27)). We conclude that also 
( s and Ce are well-behaved across the bounce. The fact that the power spectra of both (e, 
Cs and 1Z stay small across the bounce implies that perturbation theory remains valid at all 
times in our gauge, provided H(0)/Mp(0) < 1. 

We now summarize our results and draw some conclusions: 



• Having defined a non-singular bouncing cosmology, we are able to follow the behaviour 
of the various perturbation variables from beginning to end. 

• Tensor perturbations behave as expected, with h becoming constant on superhorizon 
scales. Its power spectrum can be computed by studying the associated canonical 
variable /U that turns out to be also regular throughout the bounce. The numerically 
computed power spectrum is in perfect agreement with the analytic expectations 
(including log corrections). 

• For scalar perturbations, the Bardeen potential the curvature perturbation on 
uniform dilaton hypersurfaces 1Z, and the curvature perturbation on constant density 
hypersurfaces (, all go smoothly through the bounce. 

• Both 1Z and £ stay constant on superhorizon scales, in agreement with general argu- 
ments [17], while ^ does not. 

• Provided that the ratio of the Hubble parameter to the effective Planck mass at the 
bounce (t = 0) is small, H(0)/Mp(0) < 1, 1Z and £ remain sufficiently small at all 
scales for perturbation theory to remain valid at all times in the comoving gauge. 

• By contrast, ^ becomes so large, at large scales, near the bounce that perturbation 
theory breaks down in the longitudinal gauge. 

• The canonical variable v of scalar perturbations (as well as another often discussed 
variable, u) exhibits singularities at the bounce. This is not a problem since the only 
use of v is that of giving the initial normalization of the fluctuations and the final 
Bogolubov coefficient, and v is well behaved at sufficiently early or late times. 

• We are thus able to compute numerically the scalar perturbation power spectrum. 
We find that it is very similar to that of tensor perturbations and in agreement with 
the analytic results that follow by assuming a smooth behaviour of 1Z (or C) through 
the bounce (even when the background itself has discontinuous derivatives). 
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In conclusion, at least within the class of models considered in this paper, the procedure 
used to argue that, in bouncing Universes, the spectrum of adiabatic scalar perturbations 
can be much flatter than the one of tensor perturbations, appears to be unjustified. 

We are grateful to Valerio Bozza for help during the early stages of this work. 
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